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Magnetic rtoonance imaging method 



The present invention relates to a magnetic resonance imaging method wherein 
undersampled magnetic resonance signals are acquired by a receiver antenna 
5 system having spatial sensitivity profiles according to the preamble portion of 
claim 1. The invention is further directed to a magnetic resonance apparatus aipd 
a computer program product for executing the magnetic resonance method. 

Such a magnetic resonance imaging method is usually indicated as a paralel 
10 Imaging method and known from the paper by K. Pruessmann et aL in Magn, 
Reson. Med. 4Z (1999), pages 952-962 and indicated as SENSE. Parallel imaging 
(PI) techniques exploit the spatial encoding effect inherent to Inhomogeneoiis 
receiver sensitivity. Using arrays of multiple receiver coils, parallel acquisition 
permits reducing the /c-space sampling density, thus enabling significant scan tinte 
15 reduction. PI reconstruction according to the SENSE method is typicaly 
performed by image domsiin unfolding, which enables very fest processing in the 
common case of Cartesian /c-space sampling (Cartesian SENSE). The efficien ;y 
of this approach stems from the underlying assumption that the coil sensitivities ^io 
not change significantly across each voxel. While well fulfilled at high resolution, 
20 this assumption is gradually violated at low resolution, such as in spectroscopic 
imaging (SI). As a consequence, the spatial response function (SRF) in 
conventional SENSE spectroscopic Imaging is prone to imperfections, 'whii;h 
typically result in residual aliasing artifacts. 

25 To resolve the aliasing, Cartesian SENSE reconstructs each set of aliased voxels 
separately, thus decomposing the reconstruction into a series of smaler 
independent problems, which can be solved quickly. However, with ttiis 
decomposition, the aliasing is only resolved exactly at the voxel centers. Residial 
aliasing artifacts occur when the sensitivity near a coll rises at a fester rate th an 

30 tiie decay of the sine function. These artifacts are accentuated If /r-space is 
partially zero-filled: Thus, residual artifacts can arise In high-r^olution MRI, ev3n 
wiOi 50% zerofllling in /c-space, which is often used in 3D applications. 
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This problem was previously addressed by sensitivity extrapolation (see. U. Dydak 
et.al., Magn. Reson. Med. 2001; 46: p,713^722) and by the transition to modlfiisd, 
effective sensitivity values in a recursive scheme (see X. Zhao et al.,' Proc. 
5 ISMRM 2002, Abstract No. 2393) 

It is an object of the present invention to eliminate In the above mentioned para lel 
imaging method the residual aliasing artifacts arising from high resolution images 
reconstructed with parallel imaging techniques. It is a further object of the 
10 invention to provide a magnetic resonance imaging system and a computer 
program product adapted to such a method. 

This object is solved by the method as claimed in claim 1. The further objects ere 
achieved by the magnetic resonance Imaging system according to claim 14 aid 
15 ' the computer program product according to claim 1 5. 

t 

The gist of the present invention Is that residual aliasing artifacts can >© 
eliminated by the minimum-norm least-squares reconstmcHon, which is consistent 
with the "strong" reconstruction as described in the above mentioned article of 
20 K.P, PfQssmann et. al. 



25 



Further advantages of the invention are disclosed in the dependent claims and 
the following description in which an exemplified embodiment of the invention 
described with respect to the accompanying drawings, it shows 



Fig. 1 an aliased Image and the spatial response function of the marked 
pixel in the imsige, 

Fig. 2 a minimum-norm reconstruction at different iterations of the conjugate 
gradient algorithm, 
30 Fig. 3 the panels of a zoomed portion of the reconstructions. 

Fig. 4 SRF with minimum norm and conventional reconstruction. 
Fig, 5 phantom, conventional and minimum-norm reconstruction, and ' * 
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F{g. 6 diagrammatiDally a magnetic resonance imaging system 
General principles of the invention 

Consider a general MR vnaging experiment during which iV*^ complex d^ta 
samples are acquired. Each measured signal value represents the spaflal 
integral of an object-speciTic signal density function d(r) ^ modulated ^by jan 
encoding function (r) : 

where U and denotes the complex noise component 

acquired along with the resonance signal. Noise In MRI is lypically of themiai 
original, stemming from within the object as well as from the detector circuitry. 

15 The encoding functions can reflect any modulation that the objects magnetization 
undergoes in the course of the experiment Typical means of spatial encoding sre 
plane-wave modulation by linear gradient fields, sensitivity encoding by one or 
multiple receiver coils, or RF encoding by RF pulses. The encoding functions am 
also account for undesired effects, such as phase errors due to mairi field 

20 ' inhomogeneity. Incoiporating only gradient and sensitivity encoding, for Instance, 
the encoding functions read ; 

%(r} = e^""^' ^^(r) :: [2| 

where denotes the position In fc-«pace at which the value Is tateni and 

25 denotes the complex spatial sensitivity of the receiver coil with which it is 

talcen. j 
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Reconstructing an Image from the acquired data amounts to calcuJating the.values 
of the image components, which are typically called pixels. Let the index ^- coimt 
the pixel values In the desired resulting image, with denoting their tctal 
number. Typically, the calculation of the pixel values Is assumed to be linear. 
5 reflecting the linearity of the integral in Eq. [1] with respect to rf(r) . Theniima ge 
reconstruction can be written as 

p-Fm 

where the pixel values and the measured data were assembled in the vectors 
and m , respectively, which are connected through the reconstruction matrix jF 
10 with rows and columns. 



[5] 



P 



The problem of image reconstmcHon Is thus recast as the problem of choosfng t le 
entries of the matrix F . The main goal In the choice of F Is to ensure the proper 
depletion of the object. I.e. to make each pbcel value reflect signal that origlnatss 
15 firom a corresponding small volume (voxeO »n real space. In order to track tlie 
spatial weighting of MR signal In the pixel value p^ , one can Insert Eq. [l] aito Eq, 
[3] and extract the corresponding row: 

Here the expression in brackets reflects the spatial weighting of MR signal In tlie 
20 pbcel value and is thus called the spatial response function (SRF) of the pixel: 

•: m 

Each pixel In the final Image has an Individual SRF. In order to make ^ pl^lel 
represent resonance signal from the corresponding voxel, its SRF should- Mea ly 
be non-zero Inside the voxel and zero outside the voxel. However, sucb id€al 
25 SRFs can usually not be accomplished, as they must be linear comblna^6ns of 
the encoding functions (see Eq. [5]). Therefore. SRFs In MRI are genera ly 
imperfect in that they exhibit smooth rather than sharp transitions at thelvroel 
borders and ara not exactly zero outside the voxel. In the resulting imag^ ;thB$e 
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defects are perceived as artifacts (e.g. blurring or aliasing), which Omit the Image's 
usefulness, e.g.. its diagnostic value. 

In response to this problem, H is one object of the present invention to optim ze 
5 the SRF, either globally or on a pixel by pixel basis. For this puipose, the SRF is 
preferably discretlsed by sampling the encoding functions along a finite set of 
positions In 3D space. Let this set be given by 

'>.where^^{^2,...,iyrj 

Typically this set would form a regular Cartesian grid, yet it does not need to. Tie 
to discretisation may as well be performed along any finite set of positions, e.^., 
along hexagonal or radial grids. Sampling the encoding functions at yields t|ie 
encoding matrix 

Kp^^^(rp> t7p 

which has rows and iV^ columns. Using tfiis notion, the discretised'SRF of 
15 the pixel ^ Is given by the n -th row of the SRF matrix 

More generally, each row of this matrix represents the dtscretised SRF of tlie 
corresponding voxel. Note that this representation can be taken to any. level of 

20 accuracy by enhancing the number and density of the sampling positions I 

For SRF optimisation, ttie entries of F are chosen such that the Individual .SRI -s 
approach given idealised target shapes. It depends on the application whiciji SFiF 
shape is considered ideal. In temns of overall depiction fidelity, Dirac distributioj is 
25 or boxcar Ajnctions may be regarded as optimal. For the purpose of suppressing 
far-range signal contamination, Gaussian shapes may be superior or potent al 
prior knowledge about the resonance signal distribution in the object may> favojr 
other, tailored shapes. Any target SRF shape can be accommodated within tlie 
formalism explained here. Moreover, each pixel can be assigned an Individual 
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target SRF shape. This is done by a target SRF matrix T, whose 7r-th row is tfie 
discrete representation of the SRF considered ideal for the pixei with index tv; 

The difference between the terget SRFs and those actually achieved is henpe 
5 reflected by SRF^^^ - T . For pixel-by-pixel SRF optimization, this difference 
minimized row by row. In order to do so, the magnitude of each row 
^^(Bscme'"'^ Is typically expressed i^ing a vector norm along with a geneial 
linear mapping preceding the norm operation. The SRF deviation for the pixel tp 
then 

A,=|((F£-7-)4)J I« 
Where A is a matrix representing a suitable linear mapping, the subscript ^ on tlie 
right side Indicates the ^>th row vector of the matrix pn?duct in brackets, and | ^ 

denotes the chosien vector norm. In principle, any norm can be used, including tiie 
sum, maximum, and Euclidean nomis. The Euclidean norm is preferable as 
15 leads to particulariy feasible calculations. The SRF deviation is then quantified as» 

A. ^\{{FE-T)AXl = ^[{FE-T)0{FE-Tf\^ ip. 

Where 0 = AA^ and the superscript H Indicates the Hermitean adjoint, I.e. tlie 
complex conjugate transpose. 

20 Minimizing independently for each ;r using common means of linear algeb^ 
yields the reconstruction matrix 

where the superscript indicates the Moore-Penrose pseudoinverse. 

25 Reconstruction according to Eqs. [3, 12] yields optimal SRFs. One problem thlat 
may occur In applying this fomnula is ill-condrtionlng of the matrix F. ill-condltionlng 

« i 

leads to enhancement of the noise component in Eq. [4] (second term on right 

i 

i 
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hand side of Eq. [4J), which frequently reaches unacceptable levels. Therefore, 
another object of this invention Is to control noise enhancement in addition to SRF 
quality. According to another paper by KP Pruessmann et aL (Magn Res6n M^d 
200i;46:e38-651), the noise level in the pixel value ?r is given by 



^"-^ .[13fl 
where ^ denotes the noise covariance matrix of the acquired data; The 
properties of this matrix and ttie means of determining it are desc^bed }n t^e 
above mentioned paper by KP Pruessmann et al. 



THIS 
of 



10 Using these notions, the task of ensuring proper Image quality is restated 

choosing Fsuch that both and are close to minimal for each pbcel ^: 

can, for Instance, be accomplished by minimizing ttie weighted sum of squares 
theturatemis: 

minfa Al+^Ai) .I[l4i 

15 where a and are scalar, real-valued, non-negative weighting coefficlenls. 
Performing this minimization using common means of linear algebra yields tl^e 
reconstruction matrix 

F =^ TOE^(J5E0E^ +a'¥y ; [15 

For a^O and fi^O.an equivalent formula is given by 

F = a'^T(cr'E'^'¥^E+jr^6^yE^'¥^ \ [16 

The equations (3, 15, 16] determine image reconstruction that controfe- SFIF 
quafity and noise simultaneously, individually for each pixel. To use. ithesse 
formulae, several choices must be made, which are listed and described in tlie 
following; 



1 . Format of the desired image 

The format of the desired Image must be chosen, !.e. the number of pixels :in tljie 

desired image, AT^. as well as the target shape of their corresponding i^RFs, 
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SRFj'"^^(r) . Note that by specifying the target SRF shape for the pixel fc, one 
also specifies the position of the SRFs main lobe. This position will ^picslly 
correspond to the placeinent of the pixel In the image space. The main lobe \/ill 
then shift according to the pixel position. Even if the same underlying shape is 
5 chosen for all pbcels, the target SRFs will differ by this shift. Note, how6>/^r, tr at 
the target SRFs may also vary from pixel to pixel In temis of the underlying shap e. 
Note also that the spatial arrangement of the target SRFs does not need* to l)e 
regular, e.g., such that the main lobe positions form a Cartesian grid. Each bixel's 
target SRF shape may be chosen fully independently, e.g.. sue* that main' lob >s 
10 form a hexagonal or radial pattern. ' ; 

i 

In a typical embodiment, the target SRFs would be sharply localized functiors, 
e.g. DIrac, boxcar, Gaussian, Lorentzian, or Sine shapes, vwth their centers 
anranged along a Cartesian grid. Preferably, the distance between neighlb^riiig 
15 target SRF centers is diosen to be less than the smallest features of the eridodij ig 
functions. Wrth larger distances between SRF centers, the resolution of tfi4 fiii|al 
image will be suboptimal. 

Typically, the same target SRF shape is chosen for each pixel, with indt^idiial 
20 shifts corresponding to the pixel position. This has the advantage th^ iet- 
multlpllcatlon by the matrix T, as required by Eqs. [15, 16], reduces' to a 
convolution, which can be performed very efficiently using last Fourier trartsforjm 

(FFD. ; 



Especially efficient calculations result from choosing the target SRFs to be'Diric 

shapes, with their center positions identical to the set of positions us^d 

discretisation. In this case, the matrix T Is an Identity matrix (see Eq. [9]) arid 
thus be omitted. ^ 
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2. Discretisation 

The specific form of the disjaetisation must be diosen. This is done by specifying 
the total number of sampling positions, , and the positions themselves, # 
Note that ^ese positions generally lie in 3D space. Typically, they should; cover 

5 the object to be examined witii sufficient density so that the relevant features 
the resulting SRFs are captured in the SRF matrix according to Eq. pj. It 
lm|:»retive that the discretisation be sufHa'ently fine. In the common situation 
where Fourier encoding with gradient fields is the dominant encoding mechanism 
the maximum distance b^ween neighbouring $ampnng positions should typica 

10 not exceed 7p/(2k^)^ where ^^^is the maximum value that the modulus 
the vector assumed in Eq. [2]. 



iof 
is 



Regular Cartesian grids are a good choice for discretisation. Nevertheless, t^e 

positions may form any S2unpling pattem, including hexagonal or iradlal 

15 patterns, patterns with variable density, as well as patterns that are restricted jto 
irregularly shaped volumes. 

3; SRF norm and spatial weighting 

in the described embodiment, SRF deviation is quantified by a vector riorm \n 
20 combination vinth a general linear mapping. In principle, any norm and aiiy 
mapping may be used. A preferable choice is the Euclidean norm (or 2-norn^), 
leading to the reconstruction formulae stated above. 

The additional linear mapping A can be used to introduce spatial weighting, thus 
25 expressing spatially varying priorrty in SRF optimization. Two typical rationales for 
doing so are described in the following. The first one is to exploit prior knowfledcie 
about the expected signal distribution in the object. Such knowledge p&rm ts 

malqng an estimate not just of SRF deviation but of the actual signal errpr 
incurred by SRF deviations. For instance, if the support of the object Is i^nown, A 
30 . can be made diagonal, setting to zero the diagonal entries that conBsp!c)nd to 

i : 
i ' 
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positions outside the object, while the others are set to one (see J Tsao.jlEEE 
Trans Med Imaging 2001; 20: 1178-1183.). In this fashion, one can express that 
the SRF does not matter where no MR signal is expected and thus leave! mofe 
room for noise minimisation. 

5 

Often, more than just the support of the object is known. For instance, refereni^e 
or scout images can provide an estimate of the resonance signal distributicin 
inside the object. This knowledge can be exploited by setting the diagonals <ifAio 
th@ coHBsponding signal estimates. This would usually involve interpolating and 
10 potentially regridding the prior image to the grid chosen for spatial discretisatiop 
The benefit is again a more accurate measure of signal error, thus providing 
specific additional room for noise minimization, 

A yet more advanced option Is setting 6^ equal to the expected signal 

15 coviariance in the object, in this fashion, the spatial signal correlation prejsient in 
the object can be exploited as well. However, this option requires significant 
statistical knowledge about the signal covariance, which will often net be avajlabls. 

As a second reason for spatial weighting, SRF shaping may require higher SFIF 
20 fidelity in certain target regions while penalizing SRF deviations to a smaller 
degree in other regions, irrespective of prior knowledge about the specific ;sjgral 
distribution within the object For this purpose, A should be set to a digger al 
matrix of spatial weighting coefficients, which are high where high fidelity 
desired. 

23 I 

For combining prior knowledge with additional weighting, -4^'"''"' and can be 
chosen according to the previous considerations, respectively, yielding ttne ret 
weighting . ! 

30 Finally, if neither variant of spatial weighting Is desired, A should be s^t|to 
Identity matrix. 
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Apart from Euclidean and weighted Euclidean norms, the concepts of Skf 
optimization and Joint SRF and noise optimization, as described above,: api^fy 
equally to any other norm. Non-Euclidean norms will generally enhance tie 
5' computational demands significantly. However, they can also add useful featun^ 
The maximum (Chebyshev) norm, for instance, judges the quality of art Jaclual 
SRF by its maximum deviation from the ideal. This leads to so-called equi-ripple 
deviations, where the modulus maxima have the same value along each row of 
the deviation matrix SRFat^r-T . 

4. Noise covariance matrix ^ 
In order to account accurately for noise propagation into the final image, thei noibe 
covariance matrbc T should typically be constructed according to the paper ^y 
KP Pruessmann et aL (Magn Reson Med 200l;4e:638-651), bas^d on 

15 experimental noise assessment as described therein. /Mtematively, may se 
replaced by an identity matrix, thus pretending that noise in the input data is 
mutually unconrelated and of equal strength throughout As a result, the. noise 
estimation will be ress accurate, shifting the balance between noise and SF(F 
deviaQpn towards less favourable results. Another attemative is to ensure actually 

20 uncorrelated noise by decorrelating the input data in a first step, as described in 
the same paper (Magn Reson Med 2001:46:638-651). In that case. T beeombs 
an identity matrix and can be discarded without penalty. 

6. Weighting eoefnclents cr, J5 , 

25 The weighting coefficlehte a and JS , irrtroduced in Eq. [14], express the relatii/e 
importance of minimizing noise and SRF deviation, respectively. By changing their 
ratio, the balance between the two forms of error in the final image can be shifted. 
Practically either factor can be held constant, while the ratio is adapted |via t le 
other. This remaining choice is subject to the individual imaging sitoaticn. 

30 Typically, the lesulting noise and image artifact levels should be in the same ore er 



of magnitude. Preferable choices for cc and are such that cx:W and: ^ 
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approximate the actual noise covariance In the measured data and the actial 
resonance signal covariance in the object On© problem in making these choicis 
is that the strength of the resonance signal is usually not known a priori. 
Therefore, the coefficient ratio Is typically determined empirically or fi;om 
5 preliminary reconstruction with non-optimai a and P , which provide an esfim^e 
of the signal s&ength. 

Another option is to adjust the coefficient ratio iteratively, based on the' kctijal 
noise and arti^ct levels obtained in successive reconstructions. Note that for ai iy 
10 • given coefficients, both error measures can be obtained from Eqs. [10]; 

inserting F according to Eq. [15] or [16]. In particular, with iterative adjustment of 

any desired ratio of and can be obtained on a pixel by pixel has s 
Alternatively, a certain range of coefficient ratios may be covered by a series of 
reconstructions, providing a choice of ratios between A^ and A^ for each|pix<sl 
15 A final image can then be obtained by selecting the most suited instance for each 
pbcel value. Generally, image quality changes smoothly as a Hinction of p-^l^^ 
making .the precise choice of these parameters reiatlveiy uncritical. . \ 

Further crucial choices Independent of image reconstruction concern the encoding 
ZD strategy, if gradient encoding is used, the pattern of sampled positions in /c-$pa :e 
must be chosen. Typical choices include Cartesian, hexagonal, radial, or ai spiral 
patterns. However, arbitrary fr-space patterns, including random sampling, ajre 
equally consistent with the present invention. < 

25 If sensitivity encoding with one or multiple receiver coils is used, the choice a nd 
positioning of the coils is important For RF encoding the choice of the encodi ig 
scheme is likewise essential. Any combination of these encoding mechanisms is 
consistent w'rth the presented methods. Any further encoding mechanism, as lo ig 
as essentially linear, can be accommodated by incorporating the encoding^ effect 

30 into the encoding functions introduced In Eq. [1]. \ 
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Imptementation 

As stated in Eq. [3], the general reconstruction formula is given by: 

P^Fm : [1^2 

In the described embodiment the reconstruction matrix F is calculated accord ng 
5 to Eq. [IS], or [16], These two cases, labeled A and B, are distinguished below: 

A. If the numt>er of rovra of £ is less than or equal to the number of columrfs of 
the first form of F (Eqs. [15]) is preferred as it poses the smaller inversion 
problem, leading to 

where the weighting coefTiciente c?^/? have been incorporated in the matricBs 
y and 0, respectively. Using a root decomposition (e.g. the Chlolesky 
decomposition) of , 

IS Eq. [19] can be expanded and expressed equivalently as: 

p = Te{LEY i{LE)e{LEf +iyLm j 

where denotes the identity matrix. The image vector je;^ can thus be calculated 
in 4 steps: 

A1. c^Lm 

20 A2. c' ^i{LE)0{LEf +iy c 

A3. e^eiLEfc' 

AA. P^Tc" 
The. Implementation of steps A1 to A4 is descritied below. 

1 

: ! 

25 B. If the number of rows of E is larger than the number of columns of' £, tlpe 
second form of F (Eq. [16]) Is preferred, leading to: 

p = roE:^T^*£ + o^y E^'^'^m i | [22 1 
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Where, again, the weighting coefRciente a and have been incorporated in 
the matrices Hf and 6. respectively. Using L acconJIng to Eq. poj, tti'ls can 
likewise be expanded and expressed equivalently as: 

p = T[{LEf {LE)+0-J {LEfLm • ^zJ^ 

3 The image vector p can thus be calculated in 4 steps: 

B1. c = Lm 

t 

B2. c'={LEfc 
B3. c'' = {l^LEy{LE)+e*Jc' 

10 The implementation of steps B1 to B4 is described below. 

- Steps A1 and B1 are implemented by multiplying the corresponding data by tijie 

noise decorrelation matrix This creates a new set of data, c, whose entric^s 
are characterized by mutually uncorrelated noise, while each having unitinoise 

15 level. Typically, data acquired at different points in time have no or negjigibae 

noise correlation. However, data acquired simultaneously with multiple rebeiv^r 
colls often exhibit noise correlation. In this case, the application of the matrix L 
reduces to a linear recombination of the data from the different refcelv^r 
channels (see KP Pruessmann et al., Magn Reson Med 2001;46:638-651). F 3r 

20 the sake of simplicity these lineariy combined channels will be referred: to as 

noise-decorrelated channels in the following. 

- When the data are decorrelated with respect to noise, the same 'linear 

manipulation is performed on the encoding matrix. In the most typicatj case 
25 • where noise correlation is only found among different receiver coife, tqis 
manipulation affects only the sensitivity-related factor In the encoding matrpc. 
Multiplying the matrices L and E is then practically perfonmed by creating ilinear 
recombinations of the coils' sensitivity maps. These new sensiti^^: imaps 
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correspond to virtual coils characterized by receiving mutually uncorrelatid. 
unit noise. For the sake of simplicity these sensitivity maps will be referred 
as decorrelated sensitivity maps hereafter. 



to 



of 



5 - Steps A2 and B3 Involve solving an equaUon of the form x = A* b. In case 
step A2, X, A. and b represent c[, {(^LB)0{LBy +1^ ^ and 
respectively. In case of step B3, x, A atvi b represent 
{(LBy {LE)+0*y and c' . Both steps can be solved In several ways; He 
two options are described. The first way is to explicitly detemnine the inverse 
10 A* and multiply it with b. The matrix Inversion can be accomplished, 

example, by factoring A using a singular value decomposition (SVD), and 
inverting the decomposed factors. The second way is to restate the inyers 
problem asAx^b and calculate a solution procedurally. Thfe can be dohe 
using iterative algorithnris such as the conjugate gradient (CG) method, fke 
15 second way is the prefen-ed embodiment, particularly for large matrix, sizes, 

since it does not require explicit storage of the matrix A, but only requipes the 
effect of pre-multiplication wKh ^ to be applied to Intermediate vectors y. The 
latter can be done procedurally in an efHcient manner by noting the follovying 
- Performing the multiplication (LE)y is equivalent to replicating by the 
20 number of receiver coils, multiplying each replicate with one of 

decon^elated sensitivity maps, (sampled along r^), perfbiming grldding and 
fast Fourier transfonm (FFT) for each virtual channel and ^plyng 
regridding to determine the data at eadii sampled posib'on of /^-space 
each noise-decomeiated channel. ' 
25 - Performing the multiplication (L E'fy Is equivalent to gridding A^pace d ita 

for each chanriel, applying inverse FFT for each channel, mulb'plying esich 
of the resulting images by the completx conjugate of the corresr^ond 



noise-deconelated sensitivity map, regridding along r^. and sumrjiing he 
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results over the virtual channels. Further deteiils of solving Ax.- 
proceduraily can be found in the paper by KP Pruessmann et al. (Wlajin 
Reson Med 2001 ;46:638-651). 

With the gridding operations, this implementation is applicable to j noi 
5 Cartesian sampling patterns, and it additionally allows for variable dpnsjty 

sampling, density weighting, and all types of fr-«pace shutters. 

- Steps A3 and B2 involve multiplication with . This can be perfonrwsd 
either ^pllcltly or procedurally as described sdiove for steps A2 and B3; StJip 

10 A3 additionally involves multiplication with 9. If 9 is a diagonal matrix, thjs ciin 

be perfbmied efficiently as an element-wise multiplication. 

t 

- Steps AA and B4 involve multiplication with T. If the target voxels ail haye the 
same shapes and differ only by a spatial shift, this multiplication amounte to a 

15. convolution. Then, steps AA and B4 can be implemented efficiently by f^ourier 

transibmiing c" to fc-space with FFT. multiplying the transfonned data byjthe k- 
space representation of the target voxel shape centered in the field o^ view, 
and Inverse Fourier transforming the result back to geometric space.: If tlie 
target voxels are not arranged along a Cartesian grid or the spatial 

20 disoBtisaUon is not Cartesian, the two Fourier operations require additional 

gridding. 



In all situations where gridding is mentioned, the usual considerations 
oversampling. and post-gridding apodization correction. 



or 



Piaetlcal examples of the invention 

Images (266x256. 220 mm FOV. 3 mm tiilckness) were acquired from a phantom 
in a Philips Intera 1.5T scanner (Philips Medical Systems B.V.. Best, Netherlands) 
30 using a gradient echo sequence and five surface coils placed aroiind Ihe 
phantom. Reconstruction was perfomied tteratiyely using the conjugate g(adient 
algorithm (see Fig. 2), due to its efficiency and its ability to handle non-C^jtes 

I \ 
\ i 
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sampling as well. Equation [3J is generally Ill-conditioned, due to non-orthogona itjr 
of the encoding funcUons. In the present method the Ill-conditioning ha4 besn 
overcome with TIkhonov regularization. i.e. by choosing 9 as an identity matrix 

When the full Af-space has been acquired, aliasing Is eliminated by Cartesian 
SENSE reconstruction (see Fig. 3a). However, aliasing artifacts remain noticeab e. 
If /c-space is 50% zero-filled (see Fig. 3c). The artifacts are not eliminated jy 
applying Cartesian SENSE at a lower resolution without zero-filiin^ first (see Fg 
3d), but they are eliminated completely with the reconstruction described in tliis 
invention, which will be referred to as the minimum-norm reconstrucUon (s^e Fg. 
3b). 



Figure 4 shows sample SRFs as obtained with minimum norm (a.c) aiid 
conventional (b.d) reconstruction. With the latter, the subject pbcel in the lov/er l<ifl 

15 part of the phantom suffers significant contamination from an aliasing lobe jist 
outside the object This contamination is largely eliminated wrth the mfnlmuni noim 
approach. Contamination from the pixel's close neighbourhood Is also 
considerably reduced, as evident from the line plots on the right. These plote al so 
reveal that the SRF sidelobe spacing is reduced and the main lobe Is minutely 

20 nanrowed, con-esponding to slight resolution Improvement. These effects reflect 
that minimum norm reconstruction not only accounts for sensitivity gradients, but 
actually exploits them for fevourable SRF shaping. Figure 5 shows unflltered 
reconstnjcted images along with the phantom used. With the conventioral 
method, marked artifacts occur In the center, which stem from contamination frc 

25 the object border. These and minor further artifacts are largely eliminated wth 
minimum norm reconstruction. Note the visible resolution enhancement nMar t|re 
object border where coil sensitivify gradient are the steepest. 



Rg. 6 shows diagrammatlcally a magnetic resonance imaging s^tem In which t 
30 Invention is used. : 



le 
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The magnetic resonance imaging system Includes a set of main coils 10 where 3y 
a steady, uniform magnetic field is generated. The main colls are constructed, or 
example in such a manner that they enclose a tunnel-shaped examination ^pace. 
The patient to be examined is slid on a table into this tunnei*shaped exs^m^nati on 
5 space. The magnetic resonance imaging system also Includes a number of 
gradient colls 11, 12 whereby magnetic fields exhibiting spatial variations, riotably 
In the form of temporary gradients in individual directions, are generated so as to 
be superposed on the uniform magnetic field. The gradient coils 11, 12 are 
connected to a controllable power supply unit 21. The gradient coils 11, "jZ sre 
10 energized by application of an electric cunnent by means of the power supply u ift 
21. The strength, direction and duration of the gradients are controlled by control . 
of the power supply unit. The magnetic resonance imaging system also int;lud 3s 
transmission and receiving coils 13, 15 for generating RF excitation pulses ^nd lor 
picl<ing up the magnetic resonance signals, respectively. The transmission coil 13 
15 is preferably constructed as a body coil whereby (a part of) the object ; to be 
examined can be enclosed. The body coil is usually arranged in the malgnejc 
resonance imaging system in such a manner that the patient 30 to be examined, 
being arranged in the magnetic resonance imaging system, is enclosed by the 
body coil 13, The body coil 13 acte as a transmission aerial for the transmission of 
20 the RF exdtalion pulses and RF refbcusing pulses. Preferably, the body Coil 13 
involves a spatially uniform intensity distribution of the transmitted RF pu^es. The 
receiving coils 15 are preferably surface colls 15 which are arranged on or near 
the body of the patient 30 to be examined. Such surface coils 15 have a hiiih 
sensitivity for the reception of magnetic resonance signals which Is also spatially 
25 Inhomogeneous. This means that individual surface coils 15 are mainly s6hsiti /e 
for magnetic resonance signals originating from separate directions, i.ej from 
separate parts in space of the body of the patient to be examined. The C3ii 
sensitivity profile represents the spatiai sensitivity of the set of surface coils. Tlie 
transmission coils, notably surface coils, are connected to a demodulator 24 and 
30 the received magnetic resonance signals (MS) are demodulated by means jof t le 
demodulator 24, The demodulated magnetic resonance signals (DMS) are appliijd 
to a reconstruction unit. The reconstruction unit reconstructs the ma[gne:ic 
resonance image from the demodulated magnetic resonance signals (DMS) and 
on the basis of the coil sensitivity profile of the set of surface coils. Ttie cDil 



IBTpOOl 

Empf .zeit:03/07/2003 21:39 Empf .nr.:012 P-026 



03/07/2803 20:31 +41-58-4859495 



SPIERENBURG PARTNER 



5. 



19 



sensftivity profile has been measured In advance and is stored, for ejcam 3le 
electronically, in a memory unit which is included in the reconstruction urilt^ the 
reconstruction unit derives one or more image signals from the demodulai 
magnetic resonance signals (DMS), which image signals represent one or| m 

5 possibly successive magnetic resonance images. This means that the ; signal 
levels of the image signal of such a magnetic resonance image represent 
brightness values of the relevant magnetic resonance image. The reconstt-uction 
unit 25 In practice is preferably constructed as a digital image processing* unit |25 
which Is programined so as to reconstruct the magnetic resonance image from the 

10 demodulated magnetic resonance signals and on the basis of the coil sensitivity 
profile. The digital image processing unit 25 is notably programmed so! as to 
execute the reconstruction in conformity with the present invention. The jimaje 
signal from the reconstruction unit is applied to a monitor 26 so that the nrionf or 
can display the image information of the magnetic resonance image (images). It is 

15 also possible to store the image signal in a buffer unit 27 while awaiting furt^er 
processing, for example printing in the form of a hard copy. 

In order to form a magnetic resonance image or a series of successive majgne ic 
resonance images of the patient to be examined, the body of the patient is 

20 exposed to the magnetic field prevailing in the examination space. The steaoy, 
uniform magnetic field, i.e. the main field, orients a small excess number of the 
spins in the body of the patient to be examined in the direction of the mairi field. 
This generates a (small) net macroscopic magnetization in the body. These; spins 
are, for example nuclear spins such as of the hydrogen nuclei (protorts), but 

23 electron spins may also be concerned. The magnetisation is locally influended I >y 
application of the gradient fields. For example, the gradient coils 12 apply a 
selection gradient In order to select a more or less thin slice of the bocy. 
Subsequently, the transmission coils apply the RF excitation pulse to tiie 
examination space in which the part to be imaged of the patient to be examined is 

30 situated. The RF excitation pulse excites the spins in the selected slice, i.e. the net 
magnetization then performs a precessiorial motion about the direction of th^ main 
field. During this operation those spins are excited which have a Larmor fret^ltency 
vwthin the frequency band of the RF excitation pulse in the main field. However, it 
is also very well possible to excite the spins in a part of the body which iilmu9h 
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langer man such a thin slice; for example, the spins can be excited In arthrie- 
dimensional part which extends substantially In three directions in the body. Ailer 
the RF excitation, the spins slowly return to their initial state and the macrdsco jIc 
magnetiEation returns to Its (thermal) state of equilibrium. The relaxing spins then 
emit magnetic resonance signals. Because of the application of a read-out 
gradient and a phase encoding gradient, the magnetic resonance signals havei 
plurality of frequency components which encode the spatial positions in, or 
example the selected slice. The A>-spac8 is scanned by the magnetic resonance 
signals by application of the read-out gradients and 0ie phase encoding gradien 
Accoiding to the invention, the application of notably the phase encoding 
gradients resulte in the $ub-«amF»ling of the fr-space, relative to a predetermined 
spatial resolution of the magnetic resonance image. For es^mple, a number 
lines which is too small for the predetermined resolution of the magnebc 
resonance image, for example only half the number of lines, Is scanned inithe 
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Claims 



1. A magnetic resonance Imaging method wherein spatially encoded 
undersampled magnetic resonance signals are acquired by one or; more 
5 receiver antennae and one or more images are reconstructed from tlie 

spatially encoded undersampled magnetic resonance signals^ as a result of 
optimtEiing the spatial response function (SRF), which is defined by tite 
spatial signal response from the object to be imaged. Individually for ea<:h 
pixel- 

t 

10 2. A magnetic resonance imaging method according to claim 1, wherein ttie 
spatial encoding is provided by the spatial sensitivity profiles of the rebeivsr 
antenna system, by gradient encoding, by RF encoding, by any spatial 
variation of magnetization or precession frequency, such as | fro rn 
inhomogeneity of the main magnetic field, or by any combination of ttie 

15 above encoding mechanisms. 

3. A magnetic resonance imaging method according to claim 1 or 2, wherein 
the spatial response function is optimiz^ed individually for each pixel hy 
minimizing a cost function, which determines the deviation of the spatial 
response function from a target spatial response function specified 

20 individually for each pixel i 

i 

4. A magnetic resonance imaging method according to claim 3, wherein tt^e 
cost function is determined by a norm of the deviation of the spatial response 
function after an optional linear mapping, wherein the specific choices pf 
linear mapping and norm can be diosen individually for each pixeL ! 

25 5. A magnetic resonance imaging method according to claim 4, wherein ti^e 
norm is a p-norm. ' j 

i 

6. A magnetic resonance imaging method according claim 4 or 5, wherein din 
iterative algorittim is applied for optimization of ttie cost function. 

I i 

! \ 

h 
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7. A magnetic resonance imaging method according to one of claims .1 tj 6. 
wherein the shape of the spatial response function Is optimized jointly vrith 
the noise level in the final image, and wherein optionally the cost funcjttonlfor 
each pixel is modified to additionally include the noise statistics of the pixe 

8. A magnetic resonance imaging method according to claim 7, wherein ihs 
cost fiincHon is a sum or a sum of squares of one temi measuring SRF 
deviation, and one term measuring noise statistics, with or without weighting 
in the summation. 

9. A magnetic resonance Imaging method according to claim 7, wherein tie 
estimate of the noise level is based on the noise covariance of th^ inf ut 
data. 



10, A magnetic resonance imaging method according to claims 1 to 9, where lin 
the spatial response functions are discretized with a sufficiently hiijh 
resolution to capture all significant features of target and actual SRFs. t le 

15 latter of which is determined by the resolution of the acquired data. 

11, A magnetic resonance imaging method according to one of claims 1 1to 10, 
wherein the optimization of the spatial response function Is weighted iwitt in 
the object to be imaged. 

12. A magnetic resonance imaging method according to one of claims 1 to 11, 
20 wherein target response functions that have a distinctive peak at the sir 

respective voxel centers, and arranged in any pattern, are applied ior 
optimization of the spatial response function. 

13. A magnetic resonance imaging method as in one of claims 1 to 12,- Where 
any of the choices for the linear transform or norm used In measuring SFIF 

25 deviation, the target SRFs, the pattern and density of the discretization Is 

based on an estimate of the signal distribution or covariance In the object 3r 
In reference image data, or based on anatomical or physiological side 
information. * 
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14. A magnetic resonance imaging system comprising: 

• a static main magnet having a main magnetic fieldi i 

- at least one receiver antenna, 

- means for applying a read gnadient and other gradients, 
5 - means for measuring MR signals along a predetermined trajectory 

containing a plurality of line? In fr-space 

- a receiver antenna system for acquiring undersampled MR signals^ edch 
receiver antenna position having a spatial sensithrify profle, 

- means for reconstruction the Image from the undersampled magnetic 
10 resonance signals, and 

- means for optimization of a cost function, which determines the deviafjon of 
the spatial response function from a target spatial response functi:}n 
specified Indfviduaily for each pixel and/or its corresponding noise stsiti$tics . 

15 A computer program product stored on a computer usable medium lor 
15 fonning an image by means of a magnetic resonance method comprising a 

computer readable program means for causing the computer to control me 
execution of: 

- applying a read and other gradients, | 

- measuring MR signals along a predetemiined trajectory containing 
20 plurality of lines in /r-«pace 

- acquiring undersampled MR signals from a receiver antenna system] ea^ 
receiver antenna position having a spatial sensifivity profile, ; 

reconstructing the image from the undersampled magnetic res6nan<|;e 
signals, and i 
25 - optimization of A cost function, which detennines the deviation of the'4patlal 

response function firom a target spatial response function specifi4ci 
individually for each pixel and/or its corresponding noise statistics. 
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Summary 



A novel magnetic, resonance imaging method is described, wherein underaamp 
magnetic resonance signals are acquired by a receiver antenna system ihaving 
spatial sensitivity profiles and the image being reconstructed from he 
undersampied magnetic resonance signals and the spatial sensittvity profiles. The 
reconstruction of the image is provided by an optimization of a cost function wh eh 
accounts for any of noise statistics, signal statistics, and the spatial reiporse 
function, the latter of which is defined by the spatial signal response ftxim the 
object to be imaged, separately for each individual pixel. 



(Fig. 



< • 

h 

iBTpooi I ; 

Empf.zeit:03/07/2003 21:42 " Empf .nr.:012 P.032 



1) 



03/07/2083 20: 31 +41-55-4859495 



SPIERENBURG PARThER 



33 



Fig: 1 
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Fig. 2 



a.' Cart. SENSE b. Min. Norm c. Cart. SENSE ' d. Cart:.|SE|NS& 

at full res. at half res.. 



Full acquisition 50% zero-filled. ;• 50% Z5ro-:fHled;:; 50^-z:|;4|illeU 



Fig. 3 
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